embed_wm . txt 

function y = embed_wm(x) ; 
■■]- % y = embed_wm (x) 
% 

%File Name: embed_wm . m 

%Last Modification Date: Wed Nov 1 2000 
%File Creation Date: Wed Nov 1 2000 
%Author: Jun Tian < j t ian@digimarc . com> 

%Copyright: All software, documentation, and related files in this dis 
tribution 

% are Copyright (c) 2 00 0 Digimarc Corporation 

% 

%Change History: 
% 

%Bug Report: mail to jtian@digimarc.com 
% 

*J x = haar (x) ; 

yj [m,n] = size(x)>; 

ill 

rg bits = 0; 



for i = l:m 

for j = (n/2+1) :n 

if (x(i,j) >= 0) 

bits - [bits +] ; 

else 

bits = [bits -] ; 

end 

a = abs (x (i , j ) ) + 1 ; 

a = a - 2" (floor (log (a) /log (2))) ; 

while (a > 0) 

bits = [bits (a - 2*f loor (a/2) ) ] ; 
a = a - 2 * floor (a/2) ; 

end 

end 

end 

for i = (m/2 + 1) :m 
for j = 1: (n/2) 

if (x(i,j) >= 0) 

bits - [bits +] ; 

else 

bits - [bits -] ; 
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end 

a = abs (x(i , j ) ) + 1 ; 

a = a - 2" (floor (log (a) /log (2))) ; 

while (a > 0) 

bits = [bits (a - 2*f loor (a/2) ) J 
a = a - 2 * floor (a/2) ; 

end 

end 

end 

bits - [bits SHA-l(x(l: (m/2) ,1: (n/2) ) ) ] ; 
bits = bits (2: length (bits) ) ; 
b = arithmetic_coding (bits) ; 

k = 1; 

for i - X : m 

for j = (n/2+1) :n 

q = max_capacity(x(i, j) ) ; 

a = 1 

for 1 = l:q 

a = a*2 + b (k) ; 

k = k + 1; 

end 

if (sign(x(i,j) ) >= 0) 
x(i,j) = a; 

else 

x(i ; j) - -a; 

end 

end 

end 

for i - (m/2+1) :m 
for j = 1: (n/2) 

q = max_capacity (x(i, j) ) ; 

a = 1 

for 1 = l:q 

a = a*2 + b(k) ; 
k = k + 1; 

end 

if (sign(x(i, j) ) >= 0) 
x(i,j) = a; 

else 

x(i,j) = -a; 

end 

end 

end 
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y = ihaar (x) ; 
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haar . txt 

function y = haar (x) ; 
% y = haar (x) 



%File Name: haar.m 

%Last Modification Date: Wed Nov 1 2000 
%File Creation Date: Wed Nov 1 2000 
%Author: Jun Tian < j tian@digimarc . com> 

o. 
*© 

%Copyright: All software, documentation, and related files in this dis 
tribution 

% are Copyright (c) 2000 Digimarc Corporation 

o. 

*S 

%Change History: 

% 

%Bug Report: mail to jtian@digimarc.com 

o 
*o 

fm, n] = size (x) ; 

2(1: (m/2) , :) = ceil ( (x(l:2:m, :) + X (2 :2 :m, : ) ) /2) ; 
z ( "(tn/2 + 1) :m, : ) = x(l:2:m,:) - x(2:2:m,:)/ 

y(:,l:(n/2)) = ceil ( ( z ( : , 1 : 2 : n) + z ( : , 2 : 2 :n) ) /2) ; 



Q y ( : , (n/2+1) :n) - z(: # l:2:n) - z(:,2:2:n); 
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ihaar . txt 

function y = ihaar(x); 

% y = ihaar (x) 

% 



m 



%File Name: ihaar. m 

%Last Modification Date: Wed Nov 1 2000 
%File Creation Date: Wed Nov 1 2000 
%Author: Jun Tian < j tian@digimarc . com> 
% 

%Copyright : All software, documentation, and related files in this dis 
tribution 

% are Copyright (c) 2000 Digimarc Corporation 

% 

%Change History: 
% 

%Bug Report: mail to jtian@digimarc.com 



[m / n] = size (x) ; 

z(:,l:2:n) = floor ( (2*x 1 : (n/2) ) + X ( : , (n/2+1 ) : n) ) /2) ; 

z(:,2:2:n) = f loor ( (2*x ( : , 1 : (n/2) ) - x (:, (n/2+1 ): n) ) /2 ) ; 

y(l:2:m,:) = f loor ( (2*z (1: (m/2) f : ) + z ( (m/2+1) :m, : ) ) /2) ; 

y{2:2:m / :) - f loor ( (z*x ( 1 : (m/2 ) , : ) - z ( (m/2 + 1) :m, :)} /2) ; 
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